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■ Abstract 

Formulas are derived for solutions of many-body wave scattering 
problems by small particles in the case of acoustically soft, hard, and 
impedance particles embedded in an inhomogeneous medium. The 
case of transmission (interface) boundary conditions is also studied in 
detail. The limiting case is considered, when the size a of small particles 
tends to zero while their number tends to infinity at a suitable rate. 
Equations for the limiting effective (self-consistent) field in the medium 
are derived. The theory is based on a study of integral equations and 
asymptotics of their solutions as a — ► 0. The case of wave scattering by 
many small particles embedded in an inhomogeneous medium is also 
studied. 

PACS: 02.30.Rz; 02.30.Mv; 41.20.Jb 
MSC: 35Q60;78A40; 78A45; 78A48; 
Key words: wave scattering by many small bodies; smart materials. 

1 Introduction 

There is a large literature on wave scattering by small bodies, starting from 
Rayleigh's work (1871), [H], [5], [3]. For the problem of wave scattering 
by one body an analytical solution was found only for the bodies of special 
shapes, for example, for balls and ellipsoids. If the scatterer is small then the 
scattered field can be calculated analytically for bodies of arbitrary shapes, 
sec [9j, where this theory is presented. 

The many-body wave scattering problem was discussed in the literature 
mostly numerically, if the number of scatterers is small, or under the as- 
sumption that the influence of the waves, scattered by other particles on a 
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particular particle is negligible (see [7J , where one finds a large bibliography, 
1386 entries). This corresponds to the case when the distance d between 
neighboring particles is much larger than the wavelength A, and the charac- 
teristic size a of a small body (particle) is much smaller than A. By k = ^ 
the wave number is denoted. 

This paper is a review of the author's results. The results of papers [27] 
and [3D] are presented. The presentation follows closely the above papers. 

The basic results of this paper consist of: 

i) Derivation of analytic formulas for the scattering amplitude for the 
wave scattering problem by one small {ka <C 1) body of an arbitrary shape 
under the Dirichlet, impedance, and Neumann boundary condition (acous- 
tically soft, impedance, and hard particle), and the transmission (interface) 
boundary condition; 

ii) Solution to many-body wave scattering problem by small particles, 
embedded in an inhomogeneous medium, under the assumptions a <^ d and 
a <C A, where d is the minimal distance between neighboring particles; 

iii) Derivation of the equations for the limiting effective (self-consistent) 
field in an inhomogeneous medium in which many small particles are em- 
bedded, when a — > and the number M = M (a) of the small particles tends 
to infinity at an appropriate rate; 

iv) Derivation of linear algebraic systems for solving many-body wave 
scattering problems. These systems are not obtained by a discretization of 
boundary integral equations, and they give an efficient numerical method for 
solving many-body wave scattering problems in the case of small scatterers; 

v) Theory of wave scattering by small bodies of an arbitrary shape under 
the transmission (interface) boundary condition, and a derivation of the 
equation for the effective field in the limiting medium consisting of very 
many small bodies. 

The derivations of the results, presented in this paper, are rigorous. The 
order of the error estimates as a — > is obtained. Our methods give powerful 
numerical methods for solving many-body wave scattering problems in the 
case when the scatterers are small (see [I], [2]). 

In Sections 1-4 wave scattering by small bodies under the Dirichlet, Neu- 
mann, and impedance boundary conditions is developed. In Sections 5-7 
wave scattering by small bodies under the transmission boundary condition 
is presented. In Section 8 Conclusions are briefly stated. 

Let us formulate the wave scattering problems we deal with. First, let 
us consider a one-body scattering problem. Let D\ be a bounded domain in 
R 3 with a sufficiently smooth boundary S\. The scattering problem consists 
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of finding the solution to the problem: 

(V 2 + k 2 )u = 0m D[:=R 3 \D 1 , (1) 

rn = on Si, (2) 

u = u + v, (3) 

where 

Uo = e ika ' x , aeS 2 , (4) 

S is the unit sphere in R 3 , u is the incident field, v is the scattered field 
satisfying the radiation condition 

■7 {A I I 9V M 

v r — ikv = o - , r := \x — > oo, iv := — , (5) 
\ry or 

Tu is the boundary condition (be) of one of the following types 

Tu = Tiu = u (Dirichlet be), (6) 

Tu = T2U = un — Ci^) I m Ci < 0, ( impedance be), (7) 
where £i is a constant, N is the unit normal to Si, pointing out of D\, and 

Tu = T3U = un, ( Neumann be). (8) 

The transmission be are defined in Section 5, see (jlOOp there. 

It is well known (see, e.g., [8]) that problem (H])-© has a unique solution. 
We now assume that 

a := 0.5diamL>i, ka < 1, (9) 
and look for the solution to problem (H])-© of the form 

r e ik\x-y\ 

u(x) = u {x)+ g(x,t)(Ti(t)dt, g(x,y) := — f , (10) 

where dt is the element of the surface area of Si. One can prove that the 
unique solution to the scattering problem (dJ-Q with any of the boundary 
conditions ©-(HI) can be found in the form (|10p . and the function o\ in 
equation (|10p is uniquely defined from the boundary condition ([5]). The 
scattering amplitude A(f3, a) = A(f3, a, k) is defined by the formula 

e ikr /1\ x 

v = A(B, a, k) + o - , r ->• 00, B := -. (11) 

r \r J r 
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The equations for finding a\ are: 

g(s,t)cri(t)dt = -Uq(s), 



Si 

uon - Ciuo + A<Jl „ ai ~ Ci / g{s,t)a 1 (t)dt = 0, 



(12) 
(13) 



Si 



"OAT + 



-Ao"! — 0"! 



0. 



(14) 



respectively, for conditions ([!])- 
Acr := 2 



. The operator A is defined as follows: 
d 

* <9iV. 



-g(s,t)a 1 (t)dt. 



(15) 



Equations (fT2l) - (fT4l are uniquely solvable, but there are no analytic formulas 
for their solutions for bodies of arbitrary shapes. However, if the body D\ 
is small, one can rewrite (fTUj) as 



u{x) = u (x)+g(x,0)Q 1 + / [5(x,t) - 5 (x,0)]<ri(t)dt, (16) 

J Si 



where 



Qi := / ax(t)dt, 
'Si 



(17) 



and £ Lq is the origin. 

If ka <C 1, then we prove that 



\g(x,0)Q 1 \ > 



Si 



[g(x,t) - g{x,0)]ai(t)dt 



\x\ > a. 



(18) 



Therefore, the scattered field is determined outside D\ by a single number 
Q\. This number can be obtained analytically without solving equations 
(fl~2j) - (fl~3j) . The case (fT4"|) requires a special approach by the reason discussed 
in detail later. 

Let us give the results for equations (fT2|) and (fl~3j) first. For equation 
(fl~2l) one has 



Ql= / (Ji(t)dt = -Cn (0)[l + o(l)], a^O, 
j Si 



(19) 



where C is the electric capacitance of a perfect conductor with the shape 
D\. For equation (fT3j) one has 



Qi = -C|5i|«o(0)[l + o(l)], a^O, 



(20) 
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where |Si| is the surface area of Si. The scattering amplitude for problem 
CE])-® with r = T\ (acoustically soft particle) is 

Ai03,a) = ~[l + o(l)], (21) 

since 

«o(0) = e ika - x \ x=0 = 1. 

Therefore, in this case the scattering is isotropic and of the order 0(a), 
because the capacitance C = 0(a). 

The scattering amplitude for problem (HJ)-([3J) with r = T 2 (small impedance 
particles) is : 

A 2 (a,P) = -^[l + o(l)}, (22) 

since uo(0) = 1. 

In this case the scattering is also isotropic, and of the order 0(£\Si\). 

If £ = O(l), then A 2 = 0(a 2 ), because |5i| = 0(a 2 ). If Ci = 0(±), 
k € (0, 1), then A 2 = 0(a 2 " K ). The case k = 1 was considered in 

The scattering amplitude for problem (HI)-© with r = 1^ (acoustically 
hard particles) is 



\,i i.n) --^^( ] V/ V'v)- iT "(i (2:5) 



ika-x 



Here and below summation is understood over the repeated indices, a q = 
a ■ e q , a • e q denotes the dot product of two vectors in M 3 , p, q = 1, 2, 3, {e p } 
is an orthonormal Cartesian basis of M 3 , \Di \ is the volume of D\, (3 pq is the 
magnetic polarizability tensor defined as follows ([9], p. 62): 

fipq --=7^1 [ t p a lq {t)dt, (24) 



\Di\ JSl 

o\ q is the solution to the equation 

a lq {s) = A a lq - 2N q (s), (25) 

N q (s) = N(s) ■ e g , N = N(s) is the unit outer normal to Si at the point s, 
i.e., the normal pointing out of Z?i, and is the operator A at k = 0. For 
small bodies \\A — Aq\\ = o(ka). 

If uq(x) is an arbitrary field satisfying equation ([T]), not necessarily the 
plane wave e lka ' x , then 

A 3 a) = H£i (ikPiJ^-Pp + A^o) • (26) 



5 



The above formulas are derived in Section 2. In Section 3 we develop a 
theory for many-body wave scattering problem and derive the equations for 
effective field in the medium, in which many small particles are embedded, 
as a — > 0. 

The results, presented in this paper, are based on the earlier works of 
the author f|10]-[25]). Our presentation and some of the results are novel. 
These results and methods of their derivation differ much from those in the 
homogenization theory ([1], [6]). The differences are: 

i) no periodic structure in the problems is assumed, 

ii) the operators in our problems are non-selfadjoint and have continuous 
spectrum, 

iii) the limiting medium is not homogeneous and its parameters are not 
periodic, 

iv) the technique for passing to the limit is different from one used in 
homogenization theory. 

Let us summarize the results for one-body wave scattering. 

Theorem 1.1 The scattering amplitude for the problem ([T|-(j3]) for small 
body of an arbitrary shape are given by formulas (I25p . (1261) . (I27p . for the 
boundary conditions T\, T2, T3, respectively. 



2 Derivation of the formulas for one-body wave 
scattering problems 

Let us recall the known result (see e.g., [8]) 

i} f g(x,t)a 1 (t)dt= Aai ~ ai (27) 



dN s j Sl 

concerning the limiting value of the normal derivative of single-layer po- 
tential from outside. Let x m € D m , t S S m , S m is the surface of D m , 
a = 0.5 diam£) m . 

In this Section m = 1, and x m = is the origin. 

We assume that fa « 1, ad^ 1 <^ 1, SO I X X fyi I = d^> a. Then 

„ih\x t\ pik\x Xm\ / n \ 

■7-x r = ^ -ik(x- Xm r<t-x m ) ( x + ( ka + (2g) 

4ir\x-t\ 4ir\x-x m \ V d) 

fka 2 \ 

k\x -t\ = k\x -x m \- k(x - x m )° ■ (t - x m ) + 0[ — J , (29) 
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where 

j I _ I / _ \o X ~ X fn 

U — \X X m \, \X X m ) . — -j r, 

| x x m \ 

and 

J£z!l _i + o(5). (so) 

Let us derive estimate (JT9J) . Since \t\ < a on Si, one has 

g(s,t) =g (s,t)(l + O(ka)), 

where go(s,t) = ^\ s _ t i ■ Since uo(s) is a smooth function, one has |uo(s) — 
tio(0)| = 0(a). Consequently, equation (fT2|) can be considered as an equation 
for electrostatic charge distribution a\(t) on the surface S\ of a perfect 
conductor Di, charged to the constant potential — «o(0) (up to a small term 
of the order O(ka)). It is known that the total charge Q\ = J Si (j\(t)dt of 
this conductor is equal to 

Q 1 = -Cu {0)(l + O(ka)), (31) 

where C is the electric capacitance of the perfect conductor with the shape 
Di. 

Analytic formulas for electric capacitance C of a perfect conductor of 
an arbitrary shape, which allow to calculate C with a desired accuracy, are 
derived in [9]. For example, the zeroth approximation formula is 

C (0) = T^fL. r*=|*-.|, P2) 

JSi JSi r st 

and we assume in (|32p that eo = 1, where eo is the dielectric constant of the 
homogeneous medium in which the perfect conductor is placed. Formula 
HMD is formula (HH). If u (x) = e ika ' x , then u (0) = 1, and Q x = -C(l + 
0(ka)). In this case 



which is formula ()21 j) . 

Consider now wave scattering by an impedance particle. 
Let us derive formula (f20|) . Integrate equation (fT3"|) over 5i, use the 
divergence formula 



I Si 



/ V 2 tt ^ = -/c 2 / u dx = A; 2 |Di|n (0)[l + o(l)], (33) 
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where \D\ \ = 0(a 3 ), and the formula 

-Ci / uocfa = -Ci|5i|«o(0)[l + o(l)]. (34) 
J Si 



'Si 

Futhermore | f s g(s,t)ds\ = 0(a), so 



Ci / ds g{s,t)a 1 (i)dt = 0{aQi). (35) 
J St J Si 



Therefore, the term (|35p is negligible compared with Q\ as a — > 0. Finally, 
if ka -C 1, then g(s, t) = go(s, t) (1 + ik\s — t\ + ...), and 

d d 

g{s,t) = —g (s,t)[l + O(ka)]. (36) 



dN s av ' SiV, 
Denote by A the operator 

A,a = 2jf (37) 

It is known from the potential theory that 

I A a 1 ds = -f a^dt, 2 f ^(fe = -l, i € Si. (38) 

Therefore, 

ds A(Tl ~ ai = -Qi[l + O(feo)]. (39) 



Consequently, from fromulas ([33|) - ([39|) one gets formula p2|) . 

One can see that the wave scattering by an impedance particle is isotropic, 
and the scattered field is of the order 0(£i|Si|). Since |Si| = 0(a 2 ), one 
would have 0(Ci|Si|) = 0{a 2 - K ) if Ci = O (^), k G (0, 1). 

Consider now wave scattering by an acoustically hard small particle, i.e., 
the problem with the Neumann boundary condition. 

In this case we will prove that: 

i) The scattering is anisotropic, 

ii) It is defined not by a single number, as in the previous two cases, but 
by a tensor, 

and 

iii) The order of the scattered field is 0(a 3 ) as a — > 0, for a fixed k > 0, 
i.e., the scattered field is much smaller than in the previous two cases. 
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When one integrates over Si equation (fl~3|) . one gets 



Qi= f V 2 u dx = V 2 u (0)|£>i|[l + o(l)], a^O. 

JDi 



(40) 



Thus, Qi = 0(a 3 ). Therefore, the contribution of the term e %kx °' 1 in for- 
mula d2SJ) with 

%m — will be also of the order O(a^) and should be taken 
into account, in contrast to the previous two cases. Namely, 

ufa;) = u (x) + g(x,0) [ e" ifc/3 V(t)*, = x°. (41) 

JSi \ x \ 

One has 

f e~ ikl5 t a 1 (t)dt = Q 1 -ik/3 p [ tpU^dt, (42) 
JSi J Si 

where the terms of higher order of smallness are neglected and summation 
over index p is understood. The function o~\ solves equation f j 1 4 [) : 

o\ = Ao\ + 2uqn = Aai + 2ika q NqUo(s), s € Si (43) 

if uq{x) = e %ka ' x . 

Comparing (j43]) with ()25f) , using ()24h , and taking into account that ka -C 
1, one gets 



-ifc/Sp / t p (Ti(t)dt = -ikp p \Di\Ppq(-ikaq)uo(p)[l + 0{ka)} 

J Si ' (44) 

= -k 2 \D x \p m P p a q u Q {$)[l + 0{ka)\. 

From (dDJ, ([32]) and dm one gets formula (J23J), because V 2 u = -k 2 u - 
If Uo( a; ) is an arbitrary function, satisfying equation ([1]), then ika q in 

(|43p is replaced by and —k 2 UQ = Auo, which yields formula (|26|) . 

This completes the derivation of the formulas for the solution of scalar 

wave scattering problem by one small body on the boundary of which the 

Dirichlet, or the impedance, or the Neumann boundary condition is imposed. 



3 Many-body scattering problem 

In this Section we assume that there are M = M(a) small bodies (particles) 
D m , 1 < m < M, a = 0.5maxdiamD m , ka <C 1. The distance d = d(a) 
between neighboring bodies is much larger than a, d 3> a, but we do not 
assume that d ^> A, so there may be many small particles on the distances 
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of the order of the wavelength A. This means that our medium with the 
embedded particles is not necessarily diluted. 

We assume that the small bodies are embedded in an arbitrary large 
but finite domain D, D C M 3 , so D m C D. Denote D' := M 3 \ D and 
ft := U^ =1 D mj S m := dD m , <9ft = U^ =1 S m . By N we denote a unit normal 
to <9ft, pointing out of ft, by \D m \ the volume of the body D m is denoted. 

The scattering problem consists of finding the solution to the following 
problem 

(V 2 + k 2 )u = in R 3 \ ft, (45) 

Tu = on dft, (46) 

u = u + v, (47) 

where u$ is the incident field, satisfying equation (|45p in R 3 , for example, 
uq = e tka ' x , a G S 2 , and v is the scattered field, satisfying the radiation 
condition ©. The boundary condition (f4l)j) can be of the types ©-(HJ). 
In the case of impedance boundary condition (|7|) we assume that 

un = Cm u on S m , 1 < m < M, (48) 

so the impedance may vary from one particle to another. We assume that 

Cm = %2, KG (0,1), (49) 
a K 

where x m G D m is a point in D m , and h(x), x G D, is a given function, 
which we can choose as we wish, subject to the condition Imh(x) < 0. For 
simplicity we assume that h{x) is a continuous function. 

Let us make the following assumption about the distribution of small 
particles: if A C D is an arbitrary open subset of D, then the number 
M(A) of small particles in A, assuming the impedance boundary condition, 
is: 

A/>(A) = I N(x)dx[l + o(l)l, a -> 0, (50) 

a K J A 

where N(x) > is a given function. If the Dirichlet boundary condition is 
assumed, then 

A/"d(A) = — [ N(x)dx[l + o(l)l, a^O. (51) 
The case of the Neumann boundary condition will be considered later. 
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We look for the solution to problem (|45p - (|47p with the Dirichlet bound- 
ary condition of the form 



M f 

u = u + / g(x,t)a m (t)dt, (52) 
m=l Sm 

where cr m (t) are some functions to be determined from the boundary condi- 
tion ([15]) . It is proved in [11] that problem (|45| ) -(|47| ) has a unique solution of 
the form (j52l) . For any a m (t) function (|52l) solves equation ([45]) and satisfies 
condition (|47p . The boundary condition ()46p determines cr m uniquely. How- 
ever, if M > 1, then numerical solution of the system of integral equations 
for o~ m , 1 < m < M, which one gets from the boundary condition (|46|) . is 
practically not feasible. 

To avoid this principal difficulty we prove that the solution to scattering 
problem (|45p -(|47 p is determined by M numbers 

Qm ■= / o- m (t)dt, (53) 

rather than M functions o~ m (t). 

This is possible to prove if the particles D m are small. We derive ana- 
lytical formulas for Q m as a — > 0. 

Let us define the effective (self-consistent) field u e {x) = u£ (x), acting 
on the j—th particle, by the formula 

u e (x) := u(x) — / g(x,t)aj(t)dt, \x — Xj\ ~ a. (54) 

Physically this field acts on the j—th particle and is a sum of the incident 
field and the fields acting from all other particles: 

u e {x) = u^(x) := u (x) + / g{x,t)a m (t)dt. (55) 

Let us rewrite (j55|) as follows: 

M M „ 

u e (x) = u (x) + ^2 g(x,x m )Q m + ^2 / [g(x,t) - g(x,x m )]a m (t)dt. (56) 

We want to prove that the last sum is negligible compared with the first one 
as a — > 0. To prove this, let us give some estimates. One has \t — x m \ < a, 
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\g(x,t) - g(x,x m )\ = max lo (JQ ,0 (^) j , \g(x,x m )\ = 0(l/d). 

(57) 

Therefore, if \x — Xj\ = 0(a), then 



J Sm [g(x,t) - g(x,x m )]a m (t)dt 



< 0(ad _1 + ka). (58) 



• ! Cm)Qm I 

One can also prove that 

J 1 /J 2 = 0(ka + ad- 1 ), (59) 

where J\ is the first sum in (I56D and J2 is the second sum in (]56p . Therefore, 
at any point x 6 f2' = M 3 \ one has 

M 

u e (x) = uq(x) + g(x,x m )Q m , (60) 

m=l 

where the terms of higher order of smallness are omitted. 
3.1 The case of acoustically soft particles 

If (|46p is the Dirichlet condition, then, as we have proved in Section 2 (see 
formula (j3~T1) ). one has 

Qm = ~C m U e {Xm) ■ (61) 

Thus, 

M 

u e (x) = Uq(x) - ^ g(x,x m )C m u e (x m ), x £ $7'. (62) 

771=1 

One has 

u(x) = u e (x) + o(l), a — > 0, (63) 

so the full field and effective field are practically the same. 

Let us write a linear algebraic system (LAS) for finding unknown quan- 
tities u e (x m ): 

M 

u e (xj) = u (xj) - ^ g(xj,x m )C m u e (x m ). (64) 
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If M is not very large, say M = O(10 3 ), then LAS (|64p can be solved 
numerically, and formula ([62]) can be used for calculation of u e (x). 

Consider the limiting case, when a —> 0. One can rewrite (|64p as follows: 



Ue{iq) = U Q {i q ) - ^g{iq,ip)U e {£ p ) ^ C m , (65) 

where {A p }p =1 is a union of cubes which forms a covering of D, 

max diamAp := b = b(a) ^> a, 
v 

lim b(a) = 0. (66) 

a— >0 

By Ap| we denote the volume (measure) of A p , and £ p is the center of A p , 
or a point x p in an arbitrary small body D p , located in A p . Let us assume 
that there exists the limit 

hm E ^ A ; Cm =C(C P ), £ p eA p . (67) 

For example, one may have 

C m = c{i p )a (68) 

for all m such that x m € A p , where c{x) is some function in D. If all D m 
are balls of radius a, then c(x) = 4ir. We have 

£ C m = C p aAA(A p ) = C p iV(e p )|A p |[l + o(l)], a -> 0, (69) 
so limit ([67|) exists, and 



C(C P ) = c(C p )N(Q. (70) 
From dS5J), ([5g ]l -([7D ]l one gets 

Uett q ) = u (t g )-Y,9(t q ^p)c(QN(Qu e (Z p )\A p \, l<p<P. (71) 

Linear algebraic system (I71D can be considered as the collocation method for 
solving integral equation 

u(x)=u Q (x)- g(x,y)c(y)N(y)u(y)dy. (72) 
Jd 
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It is proved in [26] that system ([71]) is uniquely solvable for all sufficiently 
small b(a), and the function 

P 

u p( x ) : = ^Zx P {x)u e {ip) (73) 
p=i 

converges in L°°(D) to the unique solution of equation (|72p . The function 
Xp(x) in ([73]) is the characteristic function of the cube A p : it is equal to 1 in 
A p and vanishes outside A p . Thus, if a — > 0, the solution to the many-body 
wave scattering problem in the case of the Dirichlet boundary condition is 
well approximated by the unique solution of the integral equation ([72]) . 

Applying the operator Lq := V 2 + k 2 to ([72]) . and using the formula 
Log(x,y) = —5(x — y), where 8(x) is the delta-function, one gets 

V 2 u + k 2 u - q(x)u = in M 3 , q(x) := c(x)N(x). (74) 

The physical conclusion is: 

If one embeds M(a) = 0(l/a) small acoustically soft particles, which are 
distributed as in (I5ip . then one creats, as a — > 0, a limiting medium, which 
is inhomogeneous, and has a refraction coefficient n 2 (x) = 1 — k~ 2 q(x). 

It is interesting from the physical point of view to note that the limit, 
as a — > 0, of the total volume of the embedded particles is zero. 

Indeed, the volume of one particle is 0(a 3 ), the total number M of the 
embedded particles is 0{o?M) = 0(a 2 ), and lim a _>.o 0(a 2 ) = 0. 

The second observation is: if (|5ip holds, then on a unit length straight 
line there are 0(^73) particles, so the distance between neighboring particles 
is d = 0(a 1 ^ 3 ). If d = 0(a 7 ) with 7 > |, then the number of the embedded 
particles in a subdomain A p is O(-h) = 0(a -37 ). In this case, for 37 > 1, 
the limit in ([69]) is C(£ p ) = lim^o c p aO{a~^) = 00. Therefore, the product 
of this limit by u remains finite only if u = in D. Physically this means 
that if the distances between neighboring perfectly soft particles are smaller 
than 0(a 1 / 3 ), namely, they are 0(a 7 ) with any 7 > ^, then u = in D. 

On the other hand, if 7 < w, then the limit C(£ p ) = 0, and u = uq in D, 
so that the embedded particles do not change, in the limit a — > 0, properties 
of the medium. 

This concludes our discussion of the scattering problem for many acous- 
tically soft particles. 
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3.2 Wave scattering by many impedance particles 

We assume now that (|l9"j) and (foT)j) hold, use the exact boundary condition 
IBJ with T = T 2 , that is, 

UeN ~ CmUe + "^"I — ~ Cm [ t)<J m {t)dt = 0, (75) 



and integrate ([75]) over 5 m in order to derive an analytical asymptotic for- 
mula for Q m = f s <r m (t)dt. 
We have 

/ u eN ds = i V 2 u e dx = 0(a 3 ), (76) 
( m u e (s)ds = h(x m )a~ K \S m \u e (x m )[l + o(l)], a -4 0, (77) 
/ i4mg ^" gm (fe = -Q w [l + (l)], a^O, (78) 



and 



Cm/ / g(s,t)<j m (t)dt = h(x m )a L R Q m = o(Q m ), < k < 1. (79) 



From (ff5]) - (j79|) one finds 

Q m = -/i(x m )a 2 - K |5 m |a- 2 u e (^ m )[l + o(l)]. (80) 

This yields the formula for the approximate solution to the wave scattering 
problem for many impedance particles: 

M 

u(x) = u (x) - a 2 ~ K ^2 g(x, x m )b m h(x m )u e (x m )[l + o(l)], (81) 

m=l 

where 

b m . — I S m I (i 

are some positive numbers which depend on the geometry of S m and are 
independent of a. For example, if all D m are balls of radius a, then b m = Air. 
A linear algebraic system for u e (x m ), analogous to (I64p . is 

M 

u e (xj) = u (xj) - a 2 ~ K ^2 g(xj,x m )bmh(x m )u e (x m ). (82) 
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The integral equation for the limiting effective field in the medium with 
embedded small particles, as a — > 0, is 



u(x) = u (x) - b I g(x,y)N(y)h(y)u(y)dy, (83) 

where 



D 



u(x) = limit e (x), (84) 

a->0 

and we have assumed in f)83[) for simplicity that b m = b for all m, that is, 
all small particles are of the same shape and size. 

Applying operator Lq = V 2 + k 2 to equation (|83p . one finds the differ- 
ential equation for the limiting effective field u(x): 

(V 2 + k 2 - bN(x)h(x))u = in M 3 , (85) 

and u satisfies condition (j47j) . 

The conclusion is: the limiting medium is inhomogeneous, and its prop- 
erties are described by the function 

q(x) := bN(x)h{x). (86) 

Since the choice of the functions N(x) > and h(x), Jmh(x) < 0, is at our 
disposal, we can create the medium with desired properties by embedding 
many small impedance particles, with suitable impedances, according to the 
distribution law (150p with a suitable N(x). The function 

1 - k- 2 q(x) = n 2 (x) (87) 

is the refraction coefficient of the limiting medium. Given a desired refrac- 
tion coefficient n 2 (x), lmn 2 (x) > 0, one can find N(x) and h(x) so that 
(|87p holds, that is, one can create a material with a desired refraction coeffi- 
cient by embedding into a given material many small particles with suitable 
boundary impedances. 

This concludes our discussion of the wave scattering problem with many 
small impedance particles. 

3.3 Wave scattering by many acoustically hard particles 

Consider now the case of acoustically hard particles, i.e., the case of Neu- 
mann boundary condition. The exact boundary integral equation for the 
function a m in this case is: 

. A m o~ m o~ m 

u eN H 7, = 0. 
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Arguing as in Section 2, see formulas ([30]) - ([31]) . one obtains 



u e (x) = u (x 



m=l 



An (r \ I ih.^ m ^ Xp ( X ™)p) du e (x m ) 
lAU e (X m ) -t- lKp m | | 



(89) 

Here we took into account that the unit vector f3 in (144p is now the vector 
i ~ . anc [ = gmjp w h ere ( x ) := x . e i s th component of 

vector x in the Euclidean orthonormal basis {e p }p =1 . 

There are three sets of unknowns in ([59"]) : u e (x m ), , and Au e (x m ), 

1 < m < M, 1 < q < 3. To obtain linear algebraic system for u e {x m ) and 
dU Q^™^ one sets x = xj in ([89]) . takes the sum in ([89]) with m ^ j. This 
yields the first set of equations for finding these unknowns. Then one takes 
derivative of equation ([89]) with respect to (x) q , sets x = Xj, and takes the 
sum in (|89p with m ^ j. This yields the second set of equations for finding 
these unknowns. Finally, one takes Laplacian of equation (|89h . sets x = Xj, 
and takes the sum in (|89p with m ^ j. This yields the third set of linear 



algebraic equations for finding u e (x m ), 



du e {x m ) 

d(x) q 



, and Au e (x r 



Passing to the limit a — > in equation (|89p. yields the equation for the 
limiting field 

u(x) = u (x)+ [ g(x,y) ( p(y)V 2 u{y) + ik ^ ^ - ^ B pq (y)) dy, (90) 



D 



dy q \x - y\ 



where p(y) and B pq (y) are defined below, see formulas ([92]) and ([93]) . 

Let us derive equation (|90p . We start by transforming the sum in (|89p. 
Let {A;}^ be a covering of Z) by cubes A;, max; diamA/ = b = b(a). We 
assume that 



b(a) » d > a, 



lim 6(a) = 0. 

a->-0 



Thus, there are many small particles D m in A;. Let xi be a point in A/. 
One has 



m=l 



^2 9{x,x m ) 
l 

L 

^2d(x,xi) 



u du e (x m ) (m) ((a) p - (g m ) p ) 

K aH Ppa i i 



Aii e (x m ) -|- iiu ■ „, 

d(x)g 



F - ^ 



(91) 
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Assume that the following limit exist: 



jjm E ^j Dm] = P(y), 02) 
a^yeAi \Ai\ 

ITa IA1 = B M> ^ 



and 



lim u e (y) = u(y), lini = ^M, lim V 2 n e (y) = V 2 n(2/)- (94) 
Then, the sum in (|9ip converges to 

/ <?(*,</) f^)V 2 nfa) + ^ y * P ~ V \ B pq {y)) dy. (95) 

Consequently, (l89j) yields in the limit a — > equation (f90l) . Equation (|90|) 
cannot be reduced to a differential equation for u(x), because ([90]) is an 
integrodifferential equation whose integrand depends on x and y. 

Let us summarize the results in the following theorem. 

Theorem 3.1. The many-body scattering problem (I45p - ()47j ) has a unique 
solution for the Dirichlet, impedance, and Neumann boundary conditions. 
The limiting effective fields in the medium obtained by embedding many small 
particles of an arbitrary shape satisfy the equations (74), (85), and (90), for 
the Dirichlet, impedance, and Neumann boundary condiions, respectively. 



4 Scattering by small particles embedded in an 
inhomogeneous medium 

Suppose that the operator V 2 + fc 2 in ([1]) and in (l4"5j) is replaced by the 
operator Lq = V 2 + k 2 n^(x), where n^(x) is a known function, 

lm.nl(x) > 0. (96) 

The function Uq(x) is the refraction coefficient of an inhomogeneous medium 
in which many small particles are embedded. The results, presented in Sec- 
tion 1-3 remain valid if one replaces function g(x, y) by the Green's function 
G(x,y), 

[V 2 + k 2 n 2 (x)]G(x,y) = -5(x-y), (97) 
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satisfying the radiation condition. We assume that 



njj(x) = 1 in D' := R 3 \ D. (98) 

The function G(x,y) is uniquely defined (see, e.g., [H]). The derivations of 
the results remain essentially the same because 

G(x,y) = g (x,y)[l + O(\x-y\)}, \x - y\ 0, (99) 

where (70(^1 y) = 4 W \ x- y \ • Estimates of G(x, y) as \x — y\ — > and as \x— y\ — > 
00 are obtained in [TT] . Smallness of particles in an inhomogeneous medium 
with refraction coefficient tlq(x) is described by the relation kn$a <C 1, where 
no := maxjieo |no(x)|, and a = maxi< m <Afdiam.D TO . 

5 Wave scattering by small bodies with transmis- 
sion (interface) boundary conditions 

There is a large literature on " homogenization" , which deals with the prop- 
erties of the medium in which other materials is distributed. Quite often it 
is assumed that the medium is periodic, and homogenization is considered 
in the framework of G-convergence ([4], [6]). In most cases, one considers 
elliptic or parabolic problems with elliptic operators positive-difinite and 
having discrete spectrum. 

A theory of wave scattering by many small particles embedded in an 
inhomogeneous medium has been developed by the author ([H]-[27]). One 
of the pratically important consequences of his theory was a derivation of 
the equation for the effective (self-consistent) field in the limiting medium, 
obtained in the limit a — Y 0, M = M(a) — > 00, where a is the characteristic 
size of a small particle, and M(a) is the total number of the embedded 
particles. 

The theory was developed in Sections 1-4 (see also papers [H]-[27]) for 
boundary conditions (be) on the surfaces of small bodies, which include the 
Dirichlet be, u\s m = 0, where S m is the surface of the m-th particle D m , 
the impedance be, Cmu\s m = u N\s m i where N is the unit normal to S m , 
pointing out of D m , £ m is the boundary impedance, and the Neumann be, 
u N \s m = 0. 

In the rest of this paper the development is presented of a similar theory 
for the transmission (interface) be: 

Pm,u^f = uTf, u + = u~ on S rn , 1 < m < M. (100) 
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Here p m is a constant, +(-) denotes the limit of J^-, from inside (outside) 
of D m . 

The physical meaning of the transmission boundary conditions is the 
continuity of the pressure and the normal component of the velocity across 
the boundaries of the discontinuity of the density. One may think about 
problem (|100j) - (|104j) (see below) as of the problem of acoustic wave scattering 
by many small bodies. 

The essential novelty of the theory, developed in this paper, is the asymp- 
totically exact, as a — > 0, treatment of the one-body and many-body scalar 
wave scattering problem in the case of small scatterers on the boundaries 
of which the transmission boundary conditions are imposed. An analytic 
explicit asymptotic formula for the field scattered by one small body is de- 
rived. An integral equation for the limiting effective field in the medium, in 
which many small bodies are embedded, is derived in the limit a —> and 
M(a) oo, where M{a) is the total number of the embedded small bodies 
(particles), and M = M(a) tends to infinity at a suitable rate as a — > 0. 

For the problem with the number M of particles not large, say, less than 
5000, our theory gives an efficient numerical method for solving many-body 
wave scattering problem. 

For the problem with M very large, say, larger than 10 5 , the solution 
to many-body wave scattering problem consists in numerical solution of the 
integral equation for the limiting field in the medium, in which small particles 
are embedded. The solution to this equation approximates the solution to 
the many-body wave scattering problem with high accuracy. 

Our approach is quite different from the approach developed in homog- 
enization theory, we do not assume periodicity in the location of the small 
scatterers. Our results are of interest also in the case when the number of 
scatterers is not large, so the homogenization theory is not applicable. 

Let us formulate the scattering problem we are treating. Below condition 
(IIOOP is assumed. Let 

M 



Let n := (J An, = R 3 \n, 



m=l 

(V 2 + k 2 )u = in ft', (101) 

(V 2 + k 2 m )u = in D m , 1 < m < M, (102) 

u = uq + v, n = e ika ' x , a G S 2 , S 2 is a unit sphere in M 3 , (103) 
fdv \ 

r (- ikvj=o(l), r — s- oo. (104) 
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We assume that p m , k and are fixed given positive constants, and the 
surfaces S m are smooth. A sufficient smoothness condition is S m € C 1 '^, [i E 
(0, 1), where S m in local coordinates is given by a continuously differentiable 
function whose first derivatives are Holder-continuous with exponent //. 
We assume that x m G D m is a point inside D m , a = idianxD m , d = 

O(as) is the distance between the neighboring particles, M(A) = ^a; m eA 1> 
is the number of particles in an arbitrary open set A, the domains D m are 
not intersecting, and 

Af(A) = — [ N(x)dx[l + o(l)}, a->0, (105) 

V J A 

where N(x) > is a function which is at our disposal, V is the volume of 
one small body, V = 0(a 3 ). If D m are balls of radius a, then V = ^f-- 

It is proved in [8] that problem (|100p - ([104p has a unique solution. 

We study wave scattering by a single small body in Section 6. In other 
words, we study in Section 6 problem f|100|) - (|104[) with M = 1, The basic 
results of this Section are formulated in Theorem 6.1. 

In section 7 wave scattering by many small bodies is considered. The 
basic results of this Section are formulated in Theorem 7.1. We always 
assume that 

ka«l, d = 0(a*). (106) 



6 Wave scattering by one small body 



Let us look for the solution to problem (|100p - (|104p with M = 1 of the form 



u(x) = u (x)+ / g(x,t)a(t)dt + x / g(x,y)u(y)dy, (107) 
Js Jd 

where S = S\, D = D±, 

ik\x-y\ 

x:=k\- k 2 , g(x, y) := — , , (108) 

47r|x — y\ 

and a(t) is to be found so that conditions (jlOOp are satisfied. For any 
a € C 0,An , \i\ € (0, 1], where C 0,Atl is the set of Holder-continuous functions 
with Holder's exponent /Ui, the solution to equation (|107p satisfies equations 
(|T0T1) and (fT02l) with M = 1, and equations (jTM]) and ([TO4]) . This is easily 
checked by a direct calculation. The second condition (jlOOp is also satisfied. 
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To satisfy the first condition in equations (jlOOp with p\ = p, one has to 
satisfy the following equation 

Aa + a Aa — a , . d , 

( p -l)uo N +p— ^— + (p - 1)^-Bu = 0, (109) 

where 

Aa = 2 J ' ^A a ( t )dt, Bu = k J g(x,y)u(y)dy, (110) 

and the well-known formulas for the limiting values of the normal derivatives 
of the single-layer potential Ta := f s g(x,t)a(t)dt on S from inside and 
outside D was used. 

In [8] one finds a proof of the following existence and uniqueness result. 
Let H 2 (D) denote the usual Sobolev space of functions twice differentiable 
in L 2 -sense. 

Proposition 1. The system of equations (|107p and (I109P for the un- 
known functions a on S and u(x) in D has a solution and this solution is 
unique in C 0,Ml x H 2 (D). 

If the solution {a, u(x)\ x& d} is found, then formula (I107p defines u = 
u(x) in M 3 . 

Let us rewrite (|109p as 

a = \A(t + 2\B 1 u + 2\u Qn , (HI) 

where 

A = — B 1 u = x—— g(x,y)u(y)dy. (112) 
l + p oN s J D 

If p G [0, oo) then A € (—1, 1). Let us now use the first assumption (|106p . 
that is, the smallness of a. One has: 

g(s,t)=g (s,t)(l + O{ka)), a -> 0; g (s,t) = 1 (113) 

47r|s — t\ 

Q e lk\s-t\ Qg^ ^ ^ ^ ^ 



<w.*,|.-4 a W ,( 1 + «^»' <^ ' < 114) 

so A = Ao(l + 0((ka) 2 )), a^0;A := A\ k=0 , (115) 

B = B (1 + O(ka)), B u = x [ g (x,y)u(y)dy, (116) 

BlU = K So ^dN^ U{y)dy{1 + °^ 2a2 )) := « fl io«(l + O(fcV)). (117) 



22 



It follows from equation (|107|) that 

e ik\x— xi\ f \ f . k \ 

u(x) = u Q (x) + -i r ( — / e" lfc/ ^'V(t)iii + —u\V\ , \x - xA » a, 

\x-xi\ \4tt Js 4-k J 

(118) 

where V\ is the volume of D = D\, V\ = vol{D\) := \D\\, u\ := u(xi), 
P := \x-xi\ ■ The point xi £ D can be chosen as we wish. For one scatterer 
it is convenient to choose the origin at the point x\ so that x\ = 0. 
We did not keep the factor e~ tk P' x in the integral over D because e~ lk @' x = 
1 + 0(ka), and 

f e- ik ^ y u(y)dy = u 1 V 1 {l + 0(ka)), a -> 0. (119) 
in 

However, it will be proved that this factor under the surface integral can not 
be dropped because 

f e- lk ^a(t)dt = [ a(t)dt-ikp p [ t p a(t)dt + 0(a A ), (120) 
Js Js Js 

where over the repeated indices here and throughout this paper summation 
is understood, and the second integral in the right-hand side of (I120p is 
0(a 3 ), as a — > 0, that is, it is of the same order of smallness as the the first 
integral Q := J s o~(t)dt. The last statement will be proved later. 
With the notations 

Q := f a(t)dt, Qa := f e~ ik ^a(t)dt, (121) 

Js Js 

the expression 

A((3,a):= C ^ + ^u 1 V 1 , V\ := V := \D\, u l :=u{x 1 ), (122) 

47T 47T 

is the scattering amplitude, a is the unit vector in the direction of the 
incident wave uq = e lka ' x , (3 is the unit vector in the direction of the scattered 
wave. 

Let us prove that 

- ikpp / t p o-(t)dt = 0(a 3 ), (123) 
Js 

and therefore, the second integral in the right-hand side of equation (|120p 

cannot be dropped. 

It follows from equation (|107p that 

u(x) ^ uq(x) + g(x,x\)Qi + xg(x,xi)u(xi)Vi, \x — xi\>d^>a, (124) 
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where ~ means asymptotic equivalence as a — >■ 0. 

Formula (|124j) can be used for calculating u(x) if two quantities Q\ and 
u\ := u{x\) are found. 

Let us derive asymptotic formulas for these quantities as a — >■ 0. Integrate 
equation (jllip over S and get 

Q = 2X u 0N ds + X / Ardt + 2A / B lU ds, (125) 

Vs Vs JS 

Use formulas fjl 13p ~ (|l 17|) . the following formula (see [9J, p. 96): 



A ads = - / ads, (126) 
and the Divergence theorem, to rewrite equation (j!25|) as 

Q = 2X V 2 u dx - XQ + 2Ax / dxV 2 x g(x,y)u(y)dy. (127) 

Since 

V 2 u = -fc 2 u ; V 2 x g(x, y) = -k 2 g(x, y) - S(x - y), (128) 
equation (|127p takes the form 

(1 + X)Q = 2AV 2 ii (xi)Vi - 2Xk 2 x / dx g(x, y)udy - 2Xx I u(x)dx. 

jd Jd Jd 

(129) 

Let us use the following estimates: 

u(x)dx = uiVi(l + o(l)), a — > 0; ui := u(x%), (L30) 



/ g(x,y)u(y)dy = / dyu(y) / dxg(x,y) = 0{a ), (131) 

D JD JD JD 

f g(x,y)dx = 0(a 2 ), Vy £ D. (132) 
Jd 

From equations ([125 ]) - ([152)1 it follows that 

Q~7^TV 1 V 2 uoi--^rV 1 u 1 , a^O, (133) 
1 + A 1 + A 

where 

V 2 u i = V 2 «oOzOU =:E1 . (134) 
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Let us now integrate equation (|107j) over D and use estimate (|130p to obtain 

u\V\ = uq\V\+ I dta(t) / g(x,t)dx + k I dyu{y) I g(x,y)dx. (135) 
Js Jd Jd Jd 

If I? is a ball of radius a, then one can easily check that 

/ g(x,t)dx ~ / go(x,t)dx = — , |t| = a, a — >• 0. (136) 
Jd Jd 3 

In general, one has 

/ g(x,y)dx = 0(a 2 ), y & D, a -> 0. (137) 



If D is a ball of radius a, then equations (|135p - (|137p imply 

J* 
> 4-7ra 3 



^i = ^oi + Q 7 „ 3 + xu\0{a A ), a — > 0. (138) 
3-" 



Consequently, 

~ «oi +0(a 2 ), a^0, (139) 

because Q = 0(a 3 ). 

Indeed, from equations (|133j) and (|139p one gets 

Q ~ Vx(l-p) [V 2 n i - xnoi] , (140) 

where we took into account that 

2A 1 " P, (141) 



1 + A 



the relation u\ ~ uq\ as a — > 0, see equation (|139p . and neglected the terms 
of higher order of smallness. It follows from equation (|140p that 

Q = 0(a 3 ). (142) 

From equations (I139j) and (I140p one obtains 

u\ ~ «oi, a — > 0. (143) 

Let us now estimate Qi- One has 

Qx= a(t)dt - ik(3 p / t p a(t)dt, (144) 
25 



up to the terms of the higher order of smallness as a — >■ 0, and summation 
is understood over the repeated indices. It turns out that the integral 

I ■= t p a(t)dt (145) 
Js 

is of the same order, namely 0(a 3 ), as Q = f s a(t)dt. 
Let us check that the integral 

J ■= J dtt p-^ J d 9(t, y)u{y)dy = 0(a 4 ) 

as a — > 0, and, therefore, can be neglected compared with I. Indeed, u = 
Id m9(t,y)dy = 0{a), and f s t p dt = 0(a 3 ). Thus, J = 0(a 4 ). 
Define the function a q , q = 1,2, 3, as the unique solution to the equation 

a q = XAa q - 2XN q . (146) 

Since A = (1 — p)/(l + p), and p > 0, one concludes that A 6 (—1, 1), and 
it is known (see, for example, [9]) that the operator A is compact in L 2 (S) 
and does not have characteristic values in the interval (—1,1). This and 
the Predholm alternative imply that equation (|146p has a solution and this 
solution is unique. 

Let us prove that J s a q (t)dt = 0(a 3 ). To do this, integrate equation (|146p 
over S, take into account formula (|126p . the relation (^4 — Ao)a q = 0(a 3 ), 
and obtain 

(1 + A) / a q (t)dt = -2X [ N q dt + 0(a 3 ) = 0(a 3 ), 
Js Js 

because j s N q dt = by the Divergence theorem. 
Define the tensor 

Ppq ■= PpqW ■= Vi' 1 [ t p a q (t)dt, p,q = l, 2, 3. (147) 
Js 

This tensor is similar to the tensor f3 pq defined in [9], p. 62, by a similar 
formula with A = 1. In this case f3 pq is the magnetic polarizability tensor of 
a superconductor D placed in a homogeneous magnetic field directed along 
the unit Cartesian coordinate vector e q (see 0, p. 62). In [9] analytic 
formulas are given for calculating j3 pq with an arbitrary accuracy. 

One may neglect the term B\u in equation (jllip (because this term is 
0(a 4 )), take into account definition (I147|) . and get 

j^ tp a(t)dt = -fi pq ^V, (148) 
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where V := V\, and summation is understood over q. 
Consequently, one can rewrite ()144p as 



Qx = (1 - p)Vx[V 2 (u (xx) - xu (xx)} + ikfa-^-ppVi, /3 := _ , 

(149) 

and ( the p— th Cartesian coordinate of the vector x. 

Formula (1118P can be written as 

= ^o(x)+g(x,xi)^(l-p)[V 2 n (xi)-x^o(xi)]+iA:/3p g ^° ^ /3 p +xu (xx) 

dices, |x — »i| >> a, and 9u q}* 

^f-\y=xi, q = 1,2,3, y = (yi, 2/2,2/3) 

Formulas (1140p .( fT42]l . (1143p are valid for small D of arbitrary shape. Let 
us formulate the results of this Section in the following theorem. 

Theorem 6.1. Assume that ka <C 1, kx, k, and p are positive constants. 
Then the scattering problem (|1Q0[) - (|104|) has a unique solution. This solution 
has the form (|107p and can be calculated by formula (I150P in the region 
\x — xx\ » a up to the terms of order 0(a 4 ) as a — > 0, where a = 0.5diamD, 
x = k\ — k 2 , Vx = volD, (3 = Z~x\\ > Ppq i- s defined in equation (|147p . and 
0(a 4 ) does not depend on x. 



Here one sums over the repeated indices, \x — xx\ >> a, and 9u q}^ 



7 Wave scattering by many small bodies 

Assume that the distribution of small bodies is given by equation (|105p . and 
that there are M = M{a) non-intersecting small bodies D m of size a. For 
simplicity we assume that D m is a ball of radius a, centered at x m . There is 
an essential novel feature in the theory, developed in this paper compared 
with the one developed in [H],[IH], [21], namely, the scattered field was 
much larger, as o — > in the above papers. For example, for the impedance 
boundary condition, = on 5, the scattered field is 0(a 2 ), and for the 
Dirichlet boundary condition, u = on S, the scattered field is 0(a). 

For the Neumann boundary condition the scattered field is 0(a 3 ). We 
have the same order of smallness of the scattered field, 0(a 3 ), for the prob- 
lem with the transmission boundary condition because Vx = 0{a?). The 
basic role in Section 3 is played by formula (|150p . We assume that the 
distance d between neighboring bodies (particles) is much larger than a, 
d » a, but there can be many small particles on the wavelength, and 
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the interaction of the scattered waves (multiple scattering) is essential and 
cannot be neglected. 

This assumption effectively means that the function N(x) in (|105j) has 
to be small, N(x) « 1. Indeed, if on a segment of unit length there are 
small particles placed at a distance d between neighboring particles, then 
there are O(jj) particles on this unit segment, and O(-h) in a unit cube C\. 
since V = 0(a 3 ), by formula {TPS]) one gets 

oW)L N{x)dx = 0{ ^ 

Therefore d » a can hold only if (J c N(x)dx)^ = 0(|) << 1. 

Let us look for the (unique) solution to problem (|100p - (jl04p with 1 < 
m < M = M(a) of the form 

M . M . 

u(x) = u (x) + ^2 / g(x,t)a rn (t)dt+^2 x m g(x,y)u(y)dy. (151) 

m= l JS m m=1 J Dm 

Keeping the main terms in this equation, as a — > 0, one gets 
u(x) = u (x) + g(x,x m ) [Qm - ik ^ ( t p a m (t)dt \ + 

m= l V \ x x m\ JSm J 

M . 

+ ^ x m g(x, x m )u e (x m )V m , Q m -= a m (t)dt, a ->■ 0, (152) 

m=l , ' S ' m 

where we have used formula (|150p for the scattered field by every small 
particle, replaced uq by the effective field u e , acting on every particle, and 
took into account that /3 := (3 m := jfEf 222 ] • By (x — x m ) p the p-th component 
of vector (x — x m ) is denoted. 

The effective (self-consisted) field u e , acting on j-th particle, is defined 

as: 

M 

u e (x) = u (x)+ ^ g(x,x m ) ((1 - p m )[V 2 u e (x m ) - x m u e (x m )} + 

M 

]{x, X m )u e (x m ) Vmi \ x x j | ~ G. 

(153) 



28 



Setting x = Xj in equation (|153p . one gets a linear algebraic system for the 
unknowns uj := u e (xj),l < j < M, and du d e } X]) . Here the p— th 

component of the vector Xj, p = 1,2,3. Differentiating (j!53j) with respect 
to 

x j,p> P ~ 1)2,3, and then setting x — Xj, one obtains a linear algebraic 
system for the AM unknowns Uj and 1 < J < M, 1 < p < 3. 

This linear algebraic system one gets if one solves by a collocation method 
the following integral equation 

u(x)=u (x)+ I g(x,y)\(l -p)(V 2 - K 2 (y) + k 2 )u(y)+ 

JD 1 

ikp pq {y, X) du{y) ( f ~ y) v + (K\y) - k 2 )u(y)} N(y)dy. (154) 
dy q \x-y\ J 

In the above equation the function f3 pq (y, A) is defined as 

^(l/,A) = hm M{Ap) , 

where y = y v G A p , and tensor = is defined in (]147p . Conver- 

gence of the collocation method was proved in [24] . 

Equation (I154p is a non-local integrodifferential equation for the limiting 
effective field in the medium in which many small bodies are embedded. 

This is a novel result. The original scattering problem (|100p -( fT04j) has 
been formulated in terms of local differential operators. 

In the derivation of equation (I154p from equation (I153P we have assumed 
that p m = p does not depend on m, took into account that x 2 ^ becomes in 
the limit K 2 (y) — k 2 , and denoted by K 2 (y) a continuous function in D 
such that K 2 (x m ) = k^. As a — > the function K 2 (y) is uniquely defined 
because the set {xm}^^ becomes dense in D as a — > 0. 

To derive equation (I154p from equation (I153p we argue as follows. Con- 
sider a partition of D into a union centered at the points y p of P non- 
intersecting cubes A p , of size b(a), b(a) » d, so that each cube contains 
many small bodies, lim a ^ob(a) = 0. Let us demonstrate the passage to the 
limit a —> in the sums in equation (j!53|) using the first sum as an example. 
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Write the first sum in (|153p as 

^ ~] ffO^j -^m) (1 Pm)[^ ^e(-^m) ^m^e (•£»«,)] 
P 

p=l Xm£A p 
P 

= X! 5 ^'^^ 1 ~ ^p)[ v2n e(y P ) - K P «e(yp)]^ r (y P )|Ap|(i + 0(1)), (155) 
P =i 

where we have used formula ()105p . took into account that diamA p — > as 
a — > 0, wrote formula f|105[) as 

V ]T l = VM(A p ) = N(y p )\A p \(l + o(l)), a -+ 0, (156) 

and used the Riemann integrability of the functions involved, which holds, 
for example, if these functions are continuous. By p p we denote the value 
p(Vp)i where p{y) is a continuous function. 

The sum in f|155[) is the Riemannian sum for the integral 

/ g{x,y)(l-p(y))[V 2 u(y)-K 2 (y)u(y) + k 2 u(y)]N(y)dy. (157) 
Jd 

Similarly one treats the other sums in (|153p and obtains in the limit a — > 
equation (|154p . 

Let us formulate the results of this Section in the following theorem. 

Theorem 7.1. Assume that conditions (j!05|) and (|106p hold. Then, as 
a — > 0, the effective field, defined by equation (|153p . has a limit u(x). The 
function u(x) solves equation f|154[) . 



8 Conclusions 

In this paper analytic formulas for the scattering amplitudes for wave scat- 
tering by a single small particle are derived for various boundary conditions: 
the Dirichlet, Neumann, impedance, and transmission ones. 

The equation for the effective field in the medium, in which many small 
particles are embedded, is derived in the limit a —> 0. The physical as- 
sumptions are such that the multiple scattering effects are not negigible, 
but essential. The derivations are rigorous. 
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On the basis of the developed theory efficient numerical methods are 
proposed for solving many-body wave scattering problems in the case of 
small scatterers. 
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